/* kreisunitrechung.do --- 
 * 
 * Filename: kreisunitrechung.do
 * Description: Eigentliche Analysen auf Makro-Basis besser hier als in kreisanalyse.do
 * Author: Kai Arzheimer
 * Created: Thu Jan 13 22:17:51 2011 (+0100)
 * Version: 1.0
 * Last-Updated: Mon Mar  7 16:40:53 2011 (+0100)
 * By: Kai Arzheimer
 *     Update #: 14
 * reads: kreisunitspluscoords.dta
 * writes: 
 * Compatibility: 11.1
 * 
 */

/* Commentary: 
 * 
 * 
 * 
 */

/* Change log:
 * 
 * 
 */

/* Code: */


version 11.1
clear all
capture log close
capture cmdlog close
log using "kreisunitrechung-log.txt",replace text
cmdlog using "kreisunitrechung-cmds-log.txt",replace

* Daten lesen

use kreisunitspluscoords.dta,clear

* Korrelation PDS 2005/2002, alte L�nder ohne Berlin
corr pdswb2002 pdswb2005 if westohneberlin

*Zuwachs PDS 2005, alte L�nder ohne Berlin
summ pdsdiff0502 if westohneberlin [aw=gewichtwb2005]

* PDS 2005 als Funktion von 2002
reg pdswb2005 pdswb2002 if westohneberlin 


* PDS 2002 als Funktion von Urbanisierung und Arbeitslossigkeit in Ost und West - krass
reg pdswb2002 ost##(c.alq c.lbdichte)

* PDS 2009 
reg pdswb2009 ost##(c.alq c.lbdichte) 


*Koordinaten in km umrechnen
gen xkm=xmean/1000
gen ykm=ymean/1000

* Distanz von Saarbr�cken aus berechnen
qui summ xkm  if strmatch(terr,"*aarbr*")
local saarx = r(mean)
qui summ ykm  if strmatch(terr,"*aarbr*")
local saary = r(mean)

gen saardist = sqrt((xkm - `saarx')^2 + (ykm - `saary')^2 )
lab var saardist "Entfernung (km) von Saarbr�cken"


*Achtung, Datensatz ist zu gro� f�r das original spatial matrix Kommando. Hier gibt es eine mataisierte Version
*http://www.stata.com/statalist/archive/2010-11/msg01212.html
* spwm2

* z.B. so
*Kurzfristig nur Westen betrachten
keep if westohneberlin

*Nochmal Lafontaine-Regression
reg pdswb2005 pdswb2002 if westohneberlin
* predict rpds2005west if westohneberlin,resid
*Box-Cox-Regresssion. Erfordert strikt postive Werte
gen shiftresid=rpds2005west + 2.21
boxcox shiftresid saardist if saardist>0,model(theta)
predict boxpred
replace boxpred = boxpred - 2.21
lab var boxpred "Erwartete Residuen"
* Korrelieren sehr hoch mit tats�chlichen Residuen
corr boxpred rpds2005

* Box-Cox-Funktion
graph twoway (line boxpred saardist) (scatter rpds2005west saardist), legend(off)

*Lowess - fast noch besser
graph twoway (lowess rpds2005west saardist) (scatter rpds2005west saardist), legend(off)


*120km Luftlinie sind etwa 90 Autominuten. Extrem willkuerlich, und die Einstellung ist relevant fuer Ergebnisse
spwm2 ,name(w) xcoord(xkm) ycoord(ykm) band(0 120) standardize eigenval(e)

*R�umliche Korrelation f�r pdswb2002, pdswb2005, Residuen

spatgsa rpds,w(w) moran
spatgsa pdswb2005,w(w) moran
spatgsa pdswb2002,w(w) moran

*


exit,clear



/*  kreisunitrechung.do ends here */

